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1.  INTRODUCTION 


The  work  documented  by  this  report  was  initiated  by  the 
need  and  desire  to  better  understand  the  performance 
characteristics  of  radar  signal  processors.  State-of- 
the-art  integrated  circuit  technology  has  made  it  feasible 
to  implement  complicated  system  structures.  One  candidate 
for  use  in  radar  systems  with  heavy  ground  clutter  and 
weather  clutter  environments  is  a  two-dimensional  (2-D) 

CFAR  processor  with  clutter  map.  Raytheon  Company  performed 
an  analysis  of  this  type  of  a  processor  under  contract  with 
the  US  Army  Missile  Command  at  Redstone  Arsenal. 1  The 
Raytheon  effort  represents  a  major  contribution  to 
establishing  an  analytical  technique  for  use  with  2-D  CFAR 
signal  processors.  A  detailed  study  of  the  Raytheon  work 
was  required  so  that  each  step  of  the  highly  mathematical 
analysis  could  be  justified.  Section  2  presents  the  results 
of  the  review  of  the  theoretical  analysis.  Section  3,  4  and 
5  document  the  performance  results  obtained  after 
modification  of  the  computer  programs.  The  main  topics 
studied  were  the  effects  that  MTI  structures,  MTI 
implementations  and  window  functions  have  on  the 
probability  of  detection.  Section  6  presents  a  list  of 
proposed  extensions  to  the  analysis.  The  appendices  provide 
further  explanation  of  particular  processor  occurrences. 

2.  REVIEW  OF  THEORETICAL  PERFORMANCE  ANALYSIS 

The  block  diagram  for  the  processor  is  shown  in  Figure 
_1 .  This  is  a  linear  system  up  to  the  point  where  1  is 
calculated.  The  input  v  is  represented  as  a  column  matrix 
of  the  complex  (i.e.,  I  and  0  channels)  input  sample 
values.  It  follows  that 


?  =  A  U 

(1) 

W  =  C  X, 

(2) 

F  =  D  W  • 

(3) 

T~.  Raytheon  Company,  Automatic  Threshold  Detector  Tech¬ 
niques  ,  Final  Report  on  Contract  No.  DAAHO1-76-C-0363  for 
US  Army  Missile  Command,  ER76-4208,  15  July  1976. 
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Figure  1.  Processor  structure 


Due  to  the  linearity,  the  output  at  F  can  be  described 
by  a  Gaussian  distribution  when  Gaussian  signals  are  the 
input  to  the  system.  Also,  the  input  can  be  separated  into 
a  sum  of  components,  viz.,  ground  clutter,  weather  clutter, 
noise,  and  target  signal.  Each  component  can  be  analyzed 
separately  using  superposition.  The  prime  objective  is  to 
determine  the  variance  at  the  FFT  output.  This  is  a 
mathematically  tractable  problem  for  the  Gaussian  signals. 
Let  M  represent  the  covariance  matrix  of  v,  the 

covariance  matrix  of  If,  and  m(^)  the  covariance  matrix  of 
W.  It  follows  that 


M(]) 


- T 

=  A  M  A 


(4) 


—  (2) - (1)  —  T 

=  c  M  r 


M 


(5) 


A  similar  transformation  could  be  used  to  find  the  co- 
variance  matrix  of  F.  However,  this  isnot  necessary  be¬ 
cause  the  variance  of  each  element  of  F  is  all  that  is  re¬ 
quired.  Consequently,  the  analysis  uses  the  FFT  algorithm, 
superposition  of  the  I  and  0  signals,  and  separation  of  the 
real  and  imaginary  parts  of  the  F^  element  of  F  to  obtain 


L 

°2Rk  =  2  mij  cik  Cjk 

i/j  =  1 

L 

a2  =  V"*  m.  .  d. ,  d.. 

Ik  Z-w  i]  lk  3k 

i,j  =  1 


where  mji  terms 


are  the  elements  of  m(^) 


(6) 


(7) 


cos 


2tt_  (j-1)  k 
L 


(8) 


=  sin 


2tt  (j-1)  k 
L 


(9) 


y 


and  the  k  subscript  on  the  variance  represents  the  kth 
frequency  cell  of  the  FFT  output.  Combining  the  I  and  Q 
channel  results  yields  the  real  part  of  Fk,  i.e.,  R(Fjc)  , 
and  the  imaginary  part  of  F^,  i.e.,  I(Fk),  to  each  be  normal 
2  2 

with  variance  +  aIk  i.e., 

R(Fk)  and  I(Fk>  e  N(0,  oj^  +  a^)  .  (10) 

This  result  holds  for  the  jth  range  bin  and  the  kth 
frequency  cell  for  either  ground  clutter,  weather  clutter 
or  noise.  A  change  in  notation  is  used  to  represent  this 
feature,  viz.,  for  noise  N 

°Njk  =  °NRk  +  °NIk  .  t11) 


Similar  results  hold  for  weather  clutter,  W,  and  ground 
clutter,  g.  Thus, 


(12) 


The  magnitude  unit  will  change  the  Gaussian  distribution  of 
Fj]<  into  an  exponential  distribution  at  Zjk,  i.e., 


p(Zjk)  =_1 _  e'Zjk/2°jk 

2 

2a 

3k 


Calculation  of  the  probability  of  false  alarm  for  a  fixed 
threshold  VTk  in  the  kth  frequency  cell  yields 

PFA  =  e_Ybk 
where 

Ybk  =  VTk/2°jk  .  U5) 

When  CFAR  or  clutter  map  techniques  are  used,  the  threshold 
is  not  fixed  but  is  a  random  variable.  It  is  possible  to 
calculate  the  expected  value  (i.e.,  average  value)  of  the 
PFA  as 


10 


11 


At  a  given  range,  the  ground  clutter  backscatter  coef¬ 
ficients  are  assumed  to  be  constant  over  the  CFAR  window. 
However,  a  Wiebull  distribution  p  (o°)  is  assumed  for  the 
range  dependency  on  these  coefficients.  Thus, 


The  developments  through  Equation  (22)  are  a  terse 
parallel  of  the  Raytheon  work.  The  following  presents  an 
explanation  of  how  these  procedures  are  implemented  in  the 
Raytheon  computer  program.  The  program  is  a  hybrid  of 
equation  oriented  calculations  and  Monte  Carlo  simulation. 
The  Wiebull  statistics  of  Equation  (22)  are  evaluated  by 
Monte  Carlo  procedures.  The  selection  of  the  clutter  map  or 
2-D  CFAR  threshold  is  not  determined  by  Monte  Carlo 
methods.  Instead,  the  average  values  of  the  thresholds  from 
each  method  are  determined  and  the  largest  average  value  is 
used  to  select  the  technique.  The  mathematical  description 
of  this  is  given  below. 
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The  clutter  map  threshold  is  given  as 

M 

\  =  6k  (23) 

2  Mcr2  m  =  1 
oo 

and  the  average  value  is 


(24) 


where  o2oo  is  the  variance  in  the  bin  of  interest  and 
frequency  cell  0.  The  range  bins  below  the  bin  of  interest 
yield  a  2-D  CFAR  threshold  of 


-N 


where  represents  the  "below"  (-1)  bins  and  the 

kcn  frequency  cell.  The  average  value  is 


(26) 


Similar  results  hold  for  the  2-D  CFAR  threshold  determined 
from  the  range  bins  above  the  bin  of  interest,  i.e.,  Cv- 
and  C*.  K 

The  program  determines 


Yfck  =  Max  (Ak,  V  *  (27) 

The  selection  process  is  actually  accomplished  as 
2 

20*  Vj*  -  Max  <0^  nj.  K2  Kj  o^>  (28) 
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Once  the  threshold  is  selected,  then  the  results  of 
Equation  (21)  are  used  to  calculate  the  probability  of 
detection . 

An  actual  radar  system  with  clutter-map  and  2-D  CFAR 
would  not  use  the  average  values  as  in  Equation  (27),  but 
would  use  the  maximum  of  the  thresholds  A^,  B^  or  Cj<. 

The  threshold  selection  method  used  by  Raytheon  saves 
simulation  time  on  the  computer,  but  appears  to  be  depend¬ 
ent  on  an  intuitive  justification.  There  appears  to  be  a 
need  to  improve  this  part  of  the  analysis  and  suggestions 
are  offered  in  Section  6. 

3 .  PERFORMANCE  DEPENDENCE  ON  MTI  STRUCTURE 

A  3-pulse  canceler  and  a  narrow-notch  MTI  were  used  in 
the  analyses  performed  by  Raytheon.  It  was  desirable  to  in¬ 
vestigate  the  performance  of  other  MTI  designs;  specifical¬ 
ly,  a  2-pulse  canceler,  a  5-  and  a  9-coefficient  Chebyshev 
design  determined  by  the  Houts  and  Burlage  procedure2»3 
and  a  26-coefficient  bimodal  filter  designed  by  Holt*  to 
reject  ground  and  weather  clutter.  It  was  necessary  to  mod¬ 
ify  the  computer  program  to  incorporate  this  capability. 

The  threshold  selection  technique  was  also  refined  to 
include  a  parameter  that  had  been  omitted  in  the  original 
program.  Subsequently,  a  number  of  verification  runs  were 
made  and  the  only  differences  between  these  and  the  orig¬ 
inal  results  were  for  the  lower  numbered  frequency  cells. 
These  results  did  not  change  the  relative  performance 
results  or  the  conclusions  presented  by  Raytheon. 


2.  R.C.  Houts  and  D.W.  Burlage,  "Maximizing  the  Usable 
Bandwidth  of  MTI  Signal  Processors,"  Institute  of 
Electrical  and  Electronics  Engineers  (IEEE) 

Transactions,  on  Aerospace  and  Electronic  Systems, 

Vol.  AES-13,  No.  1,  January  1977,  pp.  48-55. 

3.  D.W.  Burlage  and  R.C.  Houts,  Design  Techniques  for 
Improved  Bandwidth  Moving  Target  Indicator  Processors 
in  Surface  Radars,  US  Army  Missile  Command,  Redstone 
Arsenal,  Alabama,  Technical  Report  RE-75-35,  1975. 

4.  B.P.  Holt,  and  R.C.  Houts,  "Multiband  FIR  Digital  Filter 
Design  Algorithm  for  Radar  Clutter  Suppression,"  Record  of 
1978  IEEE  Int.  Conf.  on  Acoustics,  Speech  and  Signal  Pro¬ 
cessing,  April  1976. 
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New  results  obtained  for  the  2-,  3-,  5-,  and  9-coef¬ 
ficient  filters  and  a  ground  clutter  environment  revealed 
practically  the  same  performance  for  all  systems  as  shown 
in  Figures  2  through  j>.  This  equivalence  of  performance  can 
be  predicted  since  all  of  the  MTI  filters  were  approximate¬ 
ly  equal  at  rejecting  the  ground  clutter.  It  might  errone¬ 
ously  be  expected  that  the  performance  of  a  CFAR  system 
with  FFT  processing  would  be  totally  independent  of  the 
MTI.  Such  faulty  reasoning  would  argue  that  a  given  fre¬ 
quency  cell  would  maintain  the  same  signal-to-interference 
ratio  with  or  without  an  MTI.  However,  the  finite  sample 
length  of  the  input  data  train  causes  sidelobe  coupling  in¬ 
to  other  frequency  cells.  The  use  of  window  functions  can 
serve  to  reduce  this  effect,  but  for  large  ground  clutter 
components  there  will  still  be  substantial  sidelobe  levels. 
The  MTI  serves  to  reduce  the  mainlobe  levels  and  conse¬ 
quently  the  sidelobes  produced  by  windowing  the  MTI 
outputs . 

A  comparison  of  the  3-pulse  MTI  to  no  MTI  is  given  in 
Figures  6,  2  and  8  as  the  spectral  width  of  the  ground 
clutter  Ts  changed.  These  results  agree  with  intuition 
since  for  a  narrow  ground  clutter  spectrum  more  of  the 
sidelobe  power  appear  at  nulls  of  the  window  function. 

When  weather  clutter  is  added  to  the  input  of  the  2-, 
3-,  5-,  and  9-coefficient  MTI  filters,  the  performance  is 
degraded  as  shown  in  Figures  9  through  12^  for  various 
weather  conditions.  The  relative  response  of  the  MTI  filter 
over  the  bandwidth  of  the  weather  has  a  major  effect  on  the 
sidelobe  contributions  and  thereby  on  the  performance  of 
other  frequency  cells.  It  is  interesting  to  note  the 
relative  improvements  in  the  performance  curves  at  fre¬ 
quencies  between  the  ground  clutter  and  weather  clutter. 

The  use  of  the  26-coefficient  bimodal  filter  is  effec¬ 
tive  in  rejecting  the  weather, provided  the  design  para¬ 
meters  are  matched  to  the  center  frequency  and  spectral 
width  of  the  clutter.  A  satisfactory  design  is  illustrated 
in  Figure  13  where  ideal  performance  is  compared  to  pre¬ 
dicted  CFAR  performance.  Comparisons  of  the  bimodal  MTI 
and  3-pulse  MTI  are  given  in  Figures  14,  15  and  16  for  var¬ 
ious  weather  conditions.  The  improved  performance  at  fre¬ 
quencies  between  the  ground  and  weather  clutters  is  a 
result  of  reducing  the  sidelobe  contributions  in  these 
cells.  This  analysis  was  performed  for  a  64-point  FFT  and 
this  requires  86  input  pulses  as  compared  to  66  pulses  when 
the  3-pulse  MTI  is  used.  When  a  32-point  FFT  is  used  with 
the  bimodal  MTI,  then  only  57  pulses  are  required.  A 
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Figure  3.  Ideal  and  predicted  probability  of  detections  with  a  3-pulse  MTI . 
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Figure  4.  Ideal  and  predicted  probability  of  detections  with  a  5-pulse  MTI 
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Figure  5.  Ideal  and  predicted  probability  of  detections  with  a  9-pulse  MTI 
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S/N  =*  10  dB,  S/C  =  -35  dB,  18-Hz  GROUND 
CLUTTER,  40-dB  TAYLOR  WEIGHTS,  16  RANGE 
BINS  IN  CFAR  WINDOW 
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Figure  12.  Effect  of  weather  on  predicted  performances  of  a  2-,  3-, 5-,  and  9-pulse 
MTI  (f«  -  945  Hz,  V  MO  ■  2000  m2) . 
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Figure  13.  Ideal  and  predicted  probability  of  d 
MTI  (no  weather). 
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FREQUENCY  CELLS 

Effect  of  weather  on  predicted  performances  of  a  3-pulse  MTI  and 
bimodal  (26-pulse)  MTI  ( f w  *  750  Hz,  ow  *  100 Hz). 
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Figure  16.  Effect  of  weather  on  predicted  performances  of  a  3-pulse  MTI  and 
bimodal  (26-pulse)  MTI  ( f w»  945  Hz,  ow  -  183.45  Hz). 


comparison  is  given  in  Figures  17  and  lj)  for  the  3-pulse, 
64-point  FFT  and  the  bimodal ,  32-point  FFT  with  and  without 
weather. 

The  relative  performance  of  the  3-pulse  and  bimodal  MTI 
structures  does  not  depend  on  the  SNR.  This  is  shown  in 
Figures  19  and  20  for  SNRs  of  0  dB  and  10  dB. 

4.  PERFORMANCE  DEPENDENCE  ON  MTI  IMPLEMENTATION 

The  MTI  used  in  radar  signal  processors  can  be 
implemented  by  either  an  Infinite  Impulse  Response  (HR)  or 
Finite  Impulse  Response  (FIR)  filter.  Due  to  transient  and 
hardware  considerations,  the  finite  impulse  response  is  the 
configuration  more  often  used.  There  are  two  major 
alternatives  for  realizing  the  finite  impulse  response 
filter.  In  radar  terminology,  these  have  become  known  by 
names  such  as  the  moving  window  and  the  fixed  window 
methods.  The  implementation  methods  allow  a  trade-off 
between  amount  of  hardware  required  and  signal-to-noise 
obtained^.  Raytheon's  original  simulation  was  based  on 
moving  window  MTI  and,  thus,  was  modified  to  provide  the 
capability  of  analysis  using  the  fixed  window  MTI. 

In  any  radar  system,  the  number  of  pulses  that  can  be 
effectively  used  is  limited  by  several  factors,  i.e.,  scan 
time,  range  walk,  target  acceleration  and  maximum  time  of 
target  coherence.  Therefore,  meaningful  performance  results 
require  maintaining  approximately  the  same  number  of 
pulses  for  different  processor  configurations.  It  has  been 
determined  that  the  moving  window  filter  produces 
residues 

KMW  -  (NP  -  N ) +1  (29) 

while  the  fixed  window  filter  produces  Kpw  residues 

Kfw  -  INTEGER  [np/n  ]  <30) 

where  Np  is  the  total  number  of  radar  pulses  available  and 
N  is  the  number  of  filter  taps. 


T~.  N.B.  Lawrence  and  J.D.  Moore,  Tradeoff  of  Hardware  and 
Siqnal-to-Noise  Ratio  in  Moving  Targets  Indicators,  US 
Army  Missile  Research  and  Development  Command,  Reds tone 
Arsenal,  Alabama,  Technical  Report  T-78-39,  AD  A054727, 
March  1978. 
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Figure  18.  Predicted  performance  comparison  of  3-Pulse  MTI ,  64-pointFFT  and  bimodal 
(26-pulse)  MTI,  32-point  FFT  (with  weather). 
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Figure  20.  Effect  of  SNR  on  predicted  performances  of  a  3-f  5-,  and  26-pulse 
MTI  (no  weather,  SNR  *  0  dB). 


Since  it  was  previously  shown  that  the  performance  is 
almost  independent  of  the  MTI  used,  the  moving  window 
versus  fixed  window  MTI  performance  analysis  was  performed 
only  on  2-pulse  and  3-pulse  cancelers,  both  commonly  used 
in  radar  systems.  Table  1  presents  a  comparison  of  radar 
pulses  available  and  residues  produced  for  various  con¬ 
figurations  of  moving  window  and  fixed  window  processors 
used  in  this  analysis. 


TABLE  1.  COMPARISON  OF  PULSES  AND  RESIDUES  FOR  MOVING 
WINDOW  AND  FIXED  WINDOW  PROCESSORS 


2-Pulse 

Canceler 

MOVING-WINDOW  } 

1  FIXED-WINDOW  1 

kmw 

NP 

KFW 

NP 

32 

33 

8 

16 

32 

33 

16 

32 

3-Pulse 

Canceler 

32 

34 

8 

24 

32 

34 

16 

48 

This  performance  analysis  is  a  computer  simulation 
which  contains  no  considerations  for  hardware,  i  .e .  „  size, 
word  length,  cost,  speed,  etc. 

The  processor,  as  analyzed,  contains  coherent 
integration,  i.e.,,-an  FFT  which  provides  3  dB  of 
signal-to-noise  improvement  for  every  factor  of  2  increase 
in  FFT  points  processed. 

The  noise  spectral  density  at  the  output  of  a  moving 
window  MTI  has  been  shaped  by  the  MTI  response,  whereas  for 
a  fixed  window  MTI  the  noise  spectrum  is  a  constant  level. 
From  an  intuitive  approach,  it  should  be  noted  that  the 
overlap  processing  of  the  moving  window  configuration 
correlates  the  input  noise  samples,  therefore# producing 
correlated  noise  at  the  output,  whereas  the  fixed  window 
configuration  generates  one  independent  noise  output  for 
several  independent  noise  inputs.  Further  explanation  is 
given  in  Appendix  A. 
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Another  important  factor  is  the  aliasing  problem 
created  by  fixed  window  filters.  Since  the  fixed  window 
filter  produces  one  output  for  N  inputs  the  sampling  rate 
has  been  reduced  by  a  factor  of  N  which  creates  aliasing 
problems.  Further  explanation  is  given  in  Appendix  B. 

In  the  performance  curves  that  follow,  it  should  be 
noted  that  the  frequency  cell  axis  corresponds  exactly  to 
the  first  16  FFT  outputs  of  the  32-point  FFT  in  the  moving 
window  system.  Only  16  are  required  due  to  symmetry.  The  8 
and  16  FFT  outputs  in  the  fixed  window  system  have  been 
rearranged  so  that  a  comparison  is  made  between  moving 
window  and  fixed  window  system  performance  for  an  input 
signal  of  a  particular  frequency.  The  rearrangement  is 
necessitated  by  two  factors.  An  obvious  difference  is 
simply  the  fact  that  an8-  or  16-  point  FFT  does  not  have  32 
outputs.  The  second  and  more  significant  difference  is  that 
while  a  movinq  window  system  has  an  FFT  with  an  input 
sample  rate  of  fs;  a  fixed  window  system  has  an  FFT  with 
an  input  sample  rate  of  fs/N.  Thus*  the  FFT  outputs  cover 
different  frequency  ranges  for  fixed  window  and  moving 
window  systems. 

The  curves  shown  in  Figures  21  and  22  are  a  2-pulse  MTI 
in  a  no  weather  environment.  In  the  lower  frequency  cells 
the  degradation  in  performance  for  fixed  window  MTI 
processor  is  due  to  fewer  FFT  points  and  signal  attenuation 
by  the  MTI  while  the  noise  is  at  a  constant  level.  For  the 
upper  cells»the  performance  degradation  is  caused  by  the 
aliasing  of  signals  so  that  they  fall  in  the  same  frequency 
cell  as  ground  clutter  and  the  reduced  FFT  length. 

The  curves  shown  in  Figures  23  and  24  are  a  2-pulse  MTI 
in  a  weather  environment.  The  same  reasons  for  performance 
degradation  hold  as  in  Figures  21  and  22.  For  the  upper 
cells,  additional  performance  degradation  is  caused  by  the 
aliasing  of  signals  so  that  they  fall  in  the  same  frequency 
cells  as  weather  clutter. 

The  curves  in  Figures  25  and  2£  are  for  a  3-pulse  MTI 
in  a  no  weather  environment.  In  both  curves  the  overall 
fixed  window  MTI  processor  performance  is  low  due  to  fewer 
FFT  points  and  to  signal  attenuation  while  the  noise  level 
remains  constant.  In  Figure  26  the  fixed  window  shows 
better  performance  than  the  moving  window  in  a  few  of  the 
upper  cells.  This  improvement  is  due  to  enhancement  of  the 
signal  by  the  MTI  while  the  noise  remains  constant. 
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FREQUENCY  CELLS 

Figure  21.  Predicted  performance  comparison  for  moving  window  2-puise  MTI 
(32-point  FFT )  and  fixed  window  2-pulse  MTI  (8-point  FFT). 
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FREQUENCY  CELLS 

Figure  22.  Predicted  performance  comparison  for  moving  window  2-pulse  MTI 
(32-point  FFT)  and  fixed  window  2-pulse  MTI  (16-point  FFT ) . 
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Fiqure  23.  Predicted  performance  comparison  for  moving  window  2-pulse  MTI 

(32-point  FFT )  and  fixed  window  2-pulse  MTI  (8-point  FFT)  in  weather 


Figure  24.  Predicted  performance  comparison  for  moving  window  2-pulse  MTI 

(32-point  FFT )  and  fixed  window  2-pulse  MTI  (16-point  FFT)  in  weather 


Figure  25.  Predicted  performance  comparison  for  moving  window  3-pulse  MTI 
(32-point  FFT)  and  fixed  window  3-pulse  MTI  (8-point  FFT). 


S/N  -  10  dB,  S/C  -  -35  dB,  18-Hz  GROUND 
CLUTTER,  40-dB  TAYLOR  WEIGHTS,  16  RANGE 
BINS  IN  CFAR  WINDOW,  3-PULSE  MTI 
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Figure  26.  Predicted  performance  comparison  for  moving  window  3-pulse  MTI 
(32-point  FFT)  and  fixed  window  3-pulse  MTI  (16-point  FFT). 


The  curves  shown  in  Figures  27  and  28  are  for  a  3-pulse 
MTI  in  a  weather  environment.  The  same  reasons  for 
performance  degradation  hold  as  in  Figures  25  and  26. 
Additional  performance  degradation  is  caused  by  the 
aliasing  of  signals  so  that  they  fall  in  the  same  frequency 
cells  as  weather  clutter. 

In  conclusion,  the  results  have  shown  that  a  moving 
window  MTI  processor  performs  considerably  better  than  a 
fixed  window  MTI  processor.  This  is  without  regard  for 
hardware  requirements.  But  using  state-of-the-art  memories 
and  multipliers  the  additional  cost  for  a  moving  window  MTI 
would  not  be  substantial  even  for  a  large  number  of  range 
cells. 

5.  PERFORMANCE  DEPENDENCE  ON  WINDOW  FUNCTION 

Typically,  in  signal  processing  systems,  a  window 
function  is  used  before  an  FFT  processor.  The  basic  effect 
of  many  window  functions  over  rectangular  is  to  reduce 
frequency  resolution  and  signal-to-noise  ratio  while 
providing  better  detection  of  tones  in  a  multitone 
environment. 

In  addition  to  multitarget  or  multitone,  a  radar 
environment  contain'-  correlated  noise  such  as  ground 
clutter,  weather  clutter  or  chaff.  Clutter  or  chaff  is 
narrowband  interference  with  a  significant  cross-section 
which,  if  not  attenuated  by  an  MTI  or  processed  with  CFAR, 
will  result  in  generating  a  substantial  number  of  false 
alarms.  One  manner  in  which  clutter  effects  processor 
performance  is  the  sidelobe  coupling  of  clutter  into  other 
FFT  frequency  cells.  Since  window  functions  are  used  to 
change  level  and  distribution  of  sidelobes  in  FFT,  the 
following  analyses  were  performed  to  determine  performance 
dependence  of  various  processor  configurations  on  window 
functions.  There  are  four  window  functions  used  in  the  , 
analysis:  1)  rectangular  (or  no  weighting),  2)  Taylor, 

3)  Kaiser,6 7 8  and  4)  Dolph-Chebyshev . ® 


6.  E.E.  Cook  and  M.  Bernfeld,  Radar  Signals:  An 
Introduction  to  Theory  and  Application,  Academic  Press, 
New  York,  196). 

7.  W.D.  Stanley,  Digital  Signal  Processing,  Prentice-Hall, 
Virginia,  1975. 

8.  Polge,  et,  al . ,  Analog-to-Digital  Converter  Evaluation 
by  Digital  Computer  Analysis  of  Recorded  Output,  US 
Army  Missile  Command,  Redstone  Arsenal,  Alabama,  Final 
Report  on  Contract  No.  DAAH01-71-C-1181 ,  1974. 
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Figure  28.  Predicted  performance  comparison  for  moving  window  3-pulse  MTI 

(32-point  FFT)  and  fixed  window  3-pulse  tni  (16-point  FFT)  in  weather 


The  performance  curves  shown  in  Figures  29  to  36  are 
for  a  signal  processor  with  an  MTI. 

Figures  29  through  32!  show  results  for  ground  clutter 
only  environment.  Since  the  MTI  has  attenuated  the  ground 
clutter#  the  noise  dominates  the  interference  with  the 
exception  of  a  few  of  the  lower  frequency  cells;  therefore, 
no  weighting  provides  best  performance  at  high  frequency 
cells  while  degrading  performance  in  the  lower  cells.  This 
is  easily  predicted  because  the  rectangular  window  function 
has  substantial  sidelobe  coupling  in  lower  cells  and  obtains 
the  best  signal-to  noise  improvement  in  noise  dominated  upper 
cells.  It  is  also  noted  that  performance  degrades  as  side- 
lobe  levels  are  lowered  because  of  the  reduction  in  the  Signal- 
to-Noise  Ratio  caused  by  the  main  lobe  of  the  window  function. 

Figures  33  through  36^  show  results  for  a  ground  and 
weather  clutter  environment.  Several  observations  can  be 
made.  The  no  weighting  system  performance  is  drastically 
reduced  due  to  sidelobe  coupling  of  weather  clutter  into 
other  frequency  cells.  For  window  functions  with  similar 
sidelobe  levels  the  performance  is  not  appreciably  differ¬ 
ent  and  the  differences  can  be  explained  by  the  mainlobe 
and  sidelobe  distributions  of  the  window  functions.  For 
this  particular  weather  input#  the  50-dB  window  functions 
provide  the  best  performance  which  is  predictable  due  to 
sidelobe  levels  that  are  needed  to  reduce  the  clutter  below 
the  noise.  That  is,  since  the  input  signal-to-noise  (S/N) 
is  10  dB  then  the  output  S/N  is  28  dB  due  to  FFT  improve¬ 
ment.  With  an  input  clutter-to-signal  (C/S)  of  35  dB  and  an 
average  MTI  clutter  attenuation  of  -5.5  dB,  the  output 
clutter-to-noise  (C/N)  is  28  dB  plus  29.5  dB  which  equals 
57.5  dB.  This  implies  that  the  clutter  would  need  to  be  re¬ 
duced  by  57.5  dB  for  the  noise  to  be  the  dominate  interference. 
Hence,  sidelobe  level  needs  to  be  down  57.5  dB.  Since  57.5 
dB  is  closer  to  50  dB  than  70  dB,  then  the  50-dB  window 
functions  should  show  best  performance. 

The  performance  curves  shown  in  Figures  37  through  44 
are  for  a  signal  processor  without  an  MTI. 

Figures  37  through  4£  show  results  for  ground  clutter 
only  environment.  Since  there  is  no  MTI  to  attenuate  the 
ground  clutter,  the  sidelobe  coupling  of  clutter  into  ail 
frequency  cells  tends  to  reduce  performance.  In  Figure  37  a 
definite  difference  in  upper  frequency  cells  in  performance 
is  noted  for  the  various  window  functions.  The  differences 
are  easily  explained  by  the  window  function  sidelobes.  The 
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Figure  29.  Predicted  performance  comparison  for  no  weighting,  40-dB  Taylor 
40-dB  Kaiser  and  40-dB  Dolph-Chebyshev  with  2-pulse  MTI. 
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Figure  30.  Predicted  performance  comparison  for  no  weighting,  50-dB  Taylor 
50-dB  Kaiser  and  50-dB  Dolph-Chebyshev  with  2-pulse  MT1. 
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WEATHER  CLUTTER 

CENTER  FREQUENCY  =  750  Hz 
CROSS  SECTION  =  200  m2 
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Figure  34.  Predicted  performance  comparison  for  no  weighting,  50-dB  Taylor, 

50-dB  Kaiser  and  50-dB  Dolph-Chebyshev  with  2-pulse  MTI  in  weather 
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Figure  38.  Predicted  performance  comparison  for  no  weighting,  50-dB  Taylor 
50-dB  Kaiser  and  50-dB  Dolph-Chebyshev  without  an  MTI . 
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FREQUENCY  CELLS 

Figure  39.  Predicted  performance  comparison  for  no  weighting,  70-dB  Taylor 
70-dB  Kaiser  and  70-dB  Dolph-Chebyshev. 
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FREQUENCY  CELLS 

Figure  41.  Predicted  performance  comparison  for  no  weighting,  40-dB  Taylor, 
40-dB  Kaiser,  and  40-dB  Dolph-Chebyshev  without  an  MTI  in  weather 
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Figure  43.  Predicted  performance  comparison  for  no  weighting,  70-dB  Taylor,  70-dB 
Kaiser,  and  70-dB  Dolph-Chebyshev  without  an  MTI  in  weather. 
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Figure  44.  Predicted  performance  comparison  for  no  weighting,  90-dB  Taylor 
90-dB  Kaiser,  and  90-dB  Dolph-Chebyshev  without  an  MTI  in  weath< 


no  weighting  has  overall  higher  sidelobes  and  the  Taylor 
and  Dolph-Chebyshev  sidelobes  remain  basically  constant  at 
40  dB  while  the  Kaiser  sidelobes  start  at  40  dB  and  de¬ 
crease  for  each  additional  sidelobe.  Thus,  the  Kaiser  has 
lower  sidelobes  at  the  lower  cells  causing  less  sidelobe 
coupling  of  ground  clutter.  For  this  particular  ground 
clutter  input,  the  70-dB  window  functions  provide  the  best 
performance  which  is  predictable  due  to  sidelobe  levels 
that  are  needed  to  reduce  clutter  below  the  noise.  That  is, 
since  the  input  S/N  is  10  dB  then  the  output  S/N  is  28  dB 
due  to  FFT  improvement.  With  an  input  clutter-to-signal 
(C/S)  of  35  dB,  the  output  clutter-to-noise  (C/N)  is  28  dB 
plus  35  dB  which  equals  63  dB.  This  implies  that  the  clut¬ 
ter  would  need  to  be  reduced  by  63  dB  for  noise  to  dominate 
the  interference.  Hence,  sidelobe  levels  need  to  be  down  63 
dB.  Since  63  dB  is  closer  to  70  dB  than  50  dB,  then  the  70- 
dB  window  functions  should  show  best  performance. 

Figures  41  to  4_4  show  results  for  a  ground  and 
weather  clutter  environment.  As  in  the  ground  only  curves 
the  performance  differences  are  due  to  the  sidelobe 
distribution  of  the  various  window  functions.  For  this 
particular  ground  and  weather  clutter  environment,  the 
70  dB  window  functions  provide  the  best  performance  results 
which  is  predictable  due  to  sidelobe  levels  that  are  needed 
to  reduce  clutter  below  the  noise.  That  is  as  previously 
determined,  for  ground  clutter  and  for  weather  clutter  the 
clutter-to-noise  (C/N)  ratio  is  63  dB  and  57.5  dB,  re¬ 
spectively.  Since  63  dB  of  sidelobes  are  required  for  noise 
to  dominate  the  interference  and  63  dB  is  closer  to  70  dB 
than  50  dB,  then  the  70  dB  window  functions  should  provide 
best  performance. 

In  conclusion,  the  results  have  shown  that  with  or 
without  an  MTI  a  window  function  with  sidelobes  sufficient 
to  reduce  clutter  below  the  noise  provides  optimum 
performance.  This  is  independent  of  which  window  function 
is  used  but  dependent  on  knowledge  of  clutter  power.  As 
shown  in  Figures  39  and  4_3  a  system  without  an  MTI  can 
obtain  equivalent  performance  to  a  system  with  an  MTI  if 
window  function  sidelobes  eliminate  clutter.  This  does  not 
take  into  consideration  necessary  FFT  wordlengths  to 
process  unattenuated  clutter  returns.  It  should  be  further 
noted  that  if  knowledge  of  clutter  power  is  limited  then 
the  best  overall  weighting,  MTI  or  no  MTI  system,  would  be 
a  window  function  with  sidelobes  that  roll-off  quickly  as 
function  of  frequency,  e.g.,  Kaiser,  since  this  would 
reduce  sidelobe  coupling  of  clutter  into  upper  frequency 
cells . 
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6. 


PROPOSED  EXTENSIONS 


There  are  a  number  of  areas  that  should  receive 
additional  studies.  Many  parameters  are  involved  in 
specifying  the  conditions  of  a  computer  run,  and  it  is  not 
feasible  to  create  a  complete  set  of  performance  curves. 

The  following  list  of  additional  runs  should  provide  more 
insight  into  the  performance  of  the  CFAR  processor. 

1.  Vary  the  signal-to-noise  ratio,  signal-to-clutter 
ratio  and  the  weather  rate. 

0 

2.  Vary  the  center  frequency  of  the  weather.  The 
analysis  is  limited  to  low  center  frequencies  and  to  "high" 
center  frequencies. 

3.  Vary  the  spectral  widths  of  the  ground  clutter  and 
the  weather  clutter. 

4.  Change  the  number  of  FFT  points  such  that  the  total 
number  of  input  pulses  remain  approximately  constant  as  the 
MTI  type  is  varied.  The  above  runs  will  not  require  program 
modification.  A  number  of  items  that  are  of  sufficient 
interest  to  warrant  the  effort  required  for  program  change 
are  listed  below. 

5.  Add  frequency-averaging  CFAR  by  determining  a 
threshold  for  range-bin  j  and  frequency  cell  k  as 


L 


E 

m  =  1 
k 


(31) 


An  alternative  is  to  only  exclude  the  spectral  magnitudes 
of  the  greatest  and  then  use  this  threshold  for  all 
frequency  cells  in  the  range  bin. 

6.  Add  the  capability  for  non-coherent  integration. 
This  appears  to  be  a  major  change  since  a  closed  form  of 
the  performance  equations  are  not  known  when  CFAR  is  used. 
However,  the  fixed  threshold  equation  for  probability  of 
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detection  can  be  obtained  from  Equation  (3-56)  of  Meyer  and 
Mayer. 9  The  fixed  threshold  probability  of  false  alarm 
calculations  are  given  by  Equation  (2-17)  of  Meyer  and 
Mayer. *  The  program  modifications  required  for  CFAR  would 
use.,  the  referenced  results  and  the  p.d.f.  descriptions 
for  the  threshold  to  calculate  expected  values  of  the 
performance  results. 

7.  A  major  change  in  the  computer  program  concept  is 
needed  to  investigate  the  effectiveness  of  the  Raytheon 
threshold  selection  technique  as  described  by  Equations 
(27)  and  (28).  One  possibility  is  to  develop  subroutines 
for  random  selection  of  values  according  to  the  p.d.f. 
descriptions  for  the  various  CFAR/clutter-map  thresholds. 
Monte-Carlo  passes  through  these  selections  and  determining 
the  largest  value  could  then  be  used  to  calculate  the 
expected  value  of  the  performance  statistics. 


9.  D.P.  Meyer,  and  H.A.  Mayer,  Radar  Target  Detection, 
Academic  Press,  New  York,  1973. 
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The  linear  system  in  Figure  A-l  where  h(t)  -  C0  6(t)  + 
Ci 6  (t-T)  +  •••  +  CN_i  6  (t-(N-l) •  T)  and  T  -  l/fin 
is  time-invariant  only  when  f0ut  equals  fin  and  the 
coefficients  CQ,  Ci ,  ...  CN_i  are  constants.  That  is, 
the  system  is  time-invariant  for  the  moving  window 
inplementation  where  fout  equals  fin  and  is 
time-varying  for  the  fixed  window  implementation  where 
f  out  equals  fin/N. 


Figure  A-l.  Linear  system 


Since  the  moving  window  MTI  is  a  time- invariant  linear 
system,  the  power  spectral  density  at  the  output  has  been 
shown  from  linear  system  analysis  to  be 


Vf)  =  Vf)’  H(f)  2  (A.U 


where  the  power  spectral  density  of  the  input  and 

|h ( f  )p  is  the  magnitude  squared  MTI  response.  Assuming  the 
input  is  white  Gaussian  noise. the  input  power  spectral 
density  is  given  by 


*xx(f)  =  a*  fr-2* 


Then,  the  output  power  spectral  density  is 


5yy(f)  =  o2 |H(f) | 2 .  (A. 3) 


Hence,  the  moving  window  MTI  output  noise  spectral  density 
is  shaped  by  the  MTI  response. 

Far  example,  If  the  MIT  is  a  2-pulse  canceler  where 
C0  ■  1  and  Ci  *  -1  it  can  be  shown  that 


1  H(f)  I  2  =  4  .sin2  ( 

irf  > 
- 2 - 

)■ 

(A. 4) 

fin  , 

1 

Substituting  (A. 2)  and  (A. 4)  into  (A.l)  we  have 


V'f) 


4^rJ  sin2 


(A. 5) 


Since  the  fixed  window  implementation  is  a  time-varying 
linear  system,  the  output  power  spectral  density  *yy( f ) 
cannot  be  determined  by  Equation  A.l.  One  technique  for 
determining  $  yy(f)  is  to  find  the  autocorrelation 
function  fyyiflff  which  faw  a  Fourier  transform  * 
pair  with  f  yy Of ) . 

For  a  fixed  window  MTI,  the  output  is  calculated  as 

y(kN)  *  N£1  Cn*M-nl,k-l,  2,  ...y  (A.6) 

n  s  0 

where  for  simplicity  T  has  been  dropped  from  the  argument  of  y(  )  and  x(  ) . 
Then  the  discrete  autocorrelation  function  +  yy(l)  is 
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♦  yytm)  ■  0,  m  jt  AN,  t  m  1,2,... 

k _ 

♦  ttn)  -  y(*N>  •  y  (  [k  +  a]  ‘  n) 


N  -  1 

N  -  1 

■  £c„ 

x  (kN  -  n). 

^  !cj x  ([k  +  j) 

n  *  0 

j  *  0  1  ’ 

N  -  1 

N  -  1 

x  (kN  -  n)0x  ^[k  +  iJ»N-j^ 

£ 

£  c"c3 

n  =  0 

j  =  o 

(A.  7) 

N  -  1 

N  -  1 

£ 

£  c"cj 

♦jCX  +  n  "  ^ 

n  ■  0 

j  =  0 

I  Letting  x(t)  De  white  Gaussian  noise  with  a  discrete 
autocorrelation  function  of 


♦„(») 


n  =  0 


=  0  ,  n  ^  0 


(A. 8) 


where  a 2  is  the  noise  variance. 


Then  the  term  x  (kN  -  n)  •  x  ( [k  +  A]  •  n  -  3)  in 
Equation  (A. 7)  is  non-zero  only  for  kN-n  *  (k  +1)  N  -  j 
or  j  •  n  +  AN.  Since  0  <  j  >  N  -  1  then  the  solution  is  l  -  0  end 
j  ■  n. 


Hence, 


(A. 9) 
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Therefore,  the  output  discrete  autocorrelation  function 
of  a  fixed  window  MTI  is  given  by 


N  -  1 

V(*N)=  E  Cn°2  =  K°2'  *  =  0 

n  =  0 

=0  /  0 


(A. 10) 


and  the  output  power  spectral  density  is  given  by 


♦yy(f)  -  *°2  (A. 11) 


Hence,  the  fixed  window  MIT  output  noise  spectral  density 
is  constant. 

For  example,  if  the  MTI  is  a  2-pulse  canceler  where 
C0=l  and  C2*-l,  then  the  power  spectral  density  is, 

N  -  1 

*yy(n  -  E  C^°2=2-°2-  (A.  12) 

n  =  0 
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APPENDIX  B 


FIXED  WINDOW  MTI  ALIASING 


MTIs  can  be  implemented  either  as  a  moving  window  or  a 
fixed  window  filter.  A  moving  window  MTI  produces  an  output 
for  each  additional  input  after  N  inputs  have  entered  the 
filter  where  N  is  the  number  of  filter  coefficients.  The 
fixed  window  produces  one  output  for  each  block  of  N  inputs 
that  enters.  Whereas  the  moving  window  MTI  provides  more 
outputs  for  integration,  the  use  of  a  fixed  window  MTI  can 
reduce  hardware  requirements. 5 

When  a  fixed  window  MTI  is  used  with  the  FFT  system, 
there  can  be  substantial  aliasing  problems.  Figure  B-1A 
shows  the  frequency  domain  representations  for  two  targets 
with  one  below  and  the  other  above  fs/3.  Figure  B-1B  shows 
a  representation  of  the  target  velocity  ambiguity  resulting 
from  a  fixed  window  three-pulse  canceler.  The  two-sided 
clutter  spectra  shown  in  Figure  B-1C  are  translated  to 
frequencies  (as  shown  in  Figure  B-1D  and  E)  where  they  can 
interfere  with  targets.  A  moving  window  MTI  would  result 
in  the  FFT  frequency  range  being  from  0  Hz  to  the  sampling 
frequency  fs.  The  fixed  window  system  effectually  reduces 
the  sample  rate  at  the  input  to  the  FFT,  thus  for  the  three- 
pulse  canceler  the  range  becomes  0  Hz  to  one-third  of  the 
original  sampling  frequency,  i.e.,  fs/3.  Based  on  this 
reduction  in  frequency  range  and  considering  the  additional 
clutter  interference  that  is  introduced,  it  is  difficult  to 
justify  the  hardware  reduction  that  could  be  obtained  by 
the  fixed  window  implementation. 
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